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ABSTRACT 

A Reduced Navier Stokes (RNS) solution technique has been successfully combined with 
the concept of partitioned geometry and mesh generation to form a very efficient 3D RNS code 
aimed at the analysis-design engineering environment. Partitioned geometry and mesh generation 
is a pre-processor to augment existing geometry and grid generation programs which allows the 
solver to (1) recluster and existing gridfile mesh lattice, and (2) perturb an existing gridfile defi- 
nition to alter the cross-sectional shape and inlet duct centerline distribution without returning 
to the external geometry and grid generator. The present results provide a quantitative validation 
of the initial value space marching 3D RNS procedure and demonstrates accurate predictions of 
the engine face flow field, with a separation present in the inlet duct as well as when vortex gen- 
erators are installed to supress flow separation. The present results also demonstrate the ability 
of the 3D RNS procedure to analyze the flow physics associated with vortex ingestion in general 
geometry ducts such as the F/A-18 inlet. At the conditions investigated, these interactions are 
basically "inviscid like", i.e. the dominant aerodynamic characteristics have there origin in inviscid 
flow theory. 


INTRODUCTION 

Modern tactical aircraft are required to be maneuverable at subsonic, transonic, and 
supersonic speeds, without giving up good cruise performance. Consequently, proper integration 
of the engine inlet with the airframe is of paramount importance. Regarding the enhancement 
of inlet performance and operation, design for optimum airframe-inlet integration has the fol- 
lowing goals: (a) to minimize approach flow angularity with respect to inlet cowl lip, (b) to deliver 
uniform, high pressure recovery flow to the inlet face, (c) to prevent or minimize vortex, wake, 
and boundary layer ingestion by the inlet over the flight envelope, (d) to reduce FOD/hot gas 
ingestion by the inlet, and finally (e) to minimize the potential for flow field interference from 
weapon carriage/firing, landing gear deployment, tanks, pods, or other hardware. The combina- 
tion of inlet design and airframe integration must not only provide high pressure recovery to 
maintain the desired thrust levels, but also generate low flow distortion consistent with stable 
engine operation. 

Engine face flow distortion is one of the most troublesome and least understood problems 
for designers of modem inlet engine systems, (Ref. 1 and 2). One issue is that there are numerous 
sources of flow field distortion that are ingested by the inlet or generated within the inlet duct 
itself. Among these sources are (a) flow separation at the cowl lip during maneuvering flight, (b) 
flow separation on the compression surfaces due to shock-wave boundary layer interactions, (c) 
spillage of the fuselage boundary layer into the inlet duct, (d) ingestion of aircraft vortices and 
wakes emanating from upstream disturbances, and (e) secondary flow and possibly flow sepa- 
ration within the inlet duct itself. Most developing aircraft have experienced one or more of these 
types of problems, particularly at high Mach numbers and/or extreme maneuver conditions, such 
that flow distortion at the engine face exceeded allowable limits. Such compatibility problems 
were encountered in the early versions of the B70, the F-l 1 1, the F-14, the MIG-25, the Tornado 
and the Airbus A300 to name a few examples. 


The effect of inlet distortion, whether it be pressure or temperature, steady or transient, is 
that the power available is reduced along with the engine-compression-system surge margin (i.e. 
the difference between the operating line and the surge line). Aeromechanical effects such as 
rotor-blade forced response and distortion effects on flutter boundaries have received less atten- 
tion, so that a consensus on importance and state-of-the-art methodology have yet to emerge. 
Stability characteristics of current high performance turbofan engines are adversely affected by 
both spatial as well as temporal distortion. According to distortion theory, there are two kinds 
of temporal distortion: (1) transients characterized by discrete frequencies (deterministic), such 
as inlet buzz during supersonic flight, vortex or wake Ingestion, and internal flow separation, and 
(2) random Auctions (stochastic), such as inlet-induced turbulence. 

The availability of present generation of super computers has opened the possibility of 
using CFD to analyze, understand, and control some of the inlet-engine flow phenomena de- 
scribed. To this end, both Full Navier Stokes and Reduced Navier Stokes analyses are being 
developed. The overall approach is to develop fast and efficient FNS and RNS computational 
techniques, and use both approaches effectively to provide overall design-guidence information 
for future aircraft concepts. The FNS code under development for design-analysis uses the highly 
efficient multi-grid approach of Vatsa and Wedan, Ref. 3, and the design-guidence RNS code 
utilizes the general velocity-decomposition approach of Briley and McDonald, Refs. 4 and 5. The 
general approach to using Computational Fluid Dynamics to study inlet-engine compatibility is- 
sues are shown schematically in Figure 1, and is based on the of using both FNS and RNS ana- 
lyses where appropriate. The Full Navier Stokes code is aimed at analysis of the inlet-forebody 
and internal duct flow field addressing such issues as (1) inlet inflow definition over a range of 
operating conditions and angles of attack, (2) detail description of the flow field in the vicinity 
of the cowl lip, (3) identification of sources of inlet flow unsteadyness such as shock boundary 
layer interactions, and internal flow separation, (4) analysis of internal flow structure within the 
transonic flow regime, (5) specialized interactions such as possible vortex breakdown within the 
inlet duct and, (6) unsteady internal interaction of a deterministic nature. The Reduced Navier 
Stokes anaysis addresses the inlet design issues such as (1) optimum duct geometry for a given 
inflow condition, (2) boundary layer control studies using devices such as vortex generators, (3) 
cross-sectional inlet duct shaping and its effect on total pressure recovery, engine face distortion, 
and transmission of disturbances through the inlet duct, and (4) distorton management control 
methods, particularily in high AOA aircraft configurations. 

This paper represents one in a series of studies on the design issues associated with inlet- 
engine compatibility, which centers on the development of CFD tools and techniques which look 
promising within a analysis-design environment, and the application of these new analysis ap- 
proaches to understand and control some of the phenomena encountered in inlet-engine com- 
patibility. Specifically, this paper deals with three topics: (1) development and validation of 
subsonic RNS computational techniques for analysis of general geometry inlet ducts such as the 
F/A-18 configuration, (2) improvement of geometry and grid generation techniques for the 
analysis-design enviroment, (3) the aerodynamics of distortion and how the transmission of dis- 
tortion through the F/A-18 inlet is effected by the geometric properties of the duct, where the 
inlet distortion arises from a vortex ingestion. 
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THEORETICAL BACKGROUND 


The reduced form of the Navier-Stokes equations originally termed parabolized Navier- 
Stokes equations (PNS) and more recently semi-elliptic, or partially parabolic are considered here 
as an initial-value space marching method for the evaluation of subsonic compressible flow with 
strong interactions and/or separation that arise from internal vortex flows. Techniques that use 
space marching with an approximate form of the Reduced Navier-Stokes equations (RNS), 
namely initial-value methods and those that require three-dimensional global iterations, have 
been considered for a number of years to predict flows in curved ducts and turbomachinery 7 blade 
casades. Unfortunately, this terminology does not identify the revelant or significant physical 
approximations made and does not distinguish these approximations from the properties of the 
solution algorithm and the differential or difference equations. In other words, different methods 
within the same "category" will in some instances give significantly different results than what is 
reported herein. 

Here, RNS methods are identified as either initial-value methods or methods which require 
global iterations. Although some type of forward-marching alogrithm may be used in both 
initial-value and global-iterations approaches, the present distinction is intended to emphasize 
that initial-value methods provide a solution after a single forward-marching integration, while 
global iteration schemes provide a solution to the system of governing equations only upon 
convergence of the iteration process. This distinction between the two categories is important, 
since conclusions reached in this paper are based on satisfying the system of governing equations 
after a single forward marching integration. A single forward marching integration is the manner 
in which computations will be accomplished in an analysis-design engineering environment. 

General Information 


Three dimensional viscous subsonic flows in complex inlet duct geometries are investigated 
by a numerical procedure which allows solution by spatial forward marching integration, utilizing 
flow approximations from the velocity-decomposition approach of Briley and McDonald, (Refs. 
4 and 5). The goal of this approach is to achieve a level of approximations that will yield accurate 
flow predictions, while reducing the labor of solution below that of the full Navier Stokes 
equations. The governing equations for this approach have been given previously for orthogonal 
coordinates, and the approach has been applied successfully to problems whose geometries can 
be fitted conveniently with orthogonal coordinates systems, (Ref. 6). However, geometries en- 
countered in typical subsonic inlet ducts cannot be treated easily using orthogonal coordinates, 
and this lead to an extension of this approach by Levy, Briley, and McDonald, Ref. 7, to treat 
ducted geometries with nonorthogonal coordinates. The nonorthogonal capability has been val- 
idated over a wide range of inlet flow conditions by Towne, (Refs. 8 and 9), and Anderson, (Ref. 
10 ). 


Partitioned Geometry and Grid Generation 

Although the analysis itself was general, the class of duct geometries that could be ana- 
lyzed was represented by superelliptic cross-sections normal to a reference line space curve having 
continous second derivatives. The description of the superelliptic cross-sections was specified by 
polynomials defined in terms of a marching parameter, r. In the present development, the ge- 
ometry formulation of Ref. 7 is extended to cover general ducted geometries define by a externally 
generated gridfile. The geometry description within the gridfile is a "ducted" geometry which has 
a variable cross-sectional area and shape, and a centerline which is curved and possibly twisted. 
In addition, the duct described by the gridfile is considered to have a defined centerline with 
continuous second derivatives. The surface geometry is described in terms of cross-sectional 
planes which lie perpendicular to the duct centerline, and thus represent the flow area at each 
streamwise station. Since the inlet duct geometry definition has been reduced to a cross section 
specification which is placed perpendicular to a centerline space curve, then a number of grid and 
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geometry pre-processing functions can be performed prior to applying the RNS flow solver. 
These pre-processing functions are shown schematically in Figure (2) and include: (1) reclustering 
the existing gridfile mesh points for more accurate solutions in regions of high shear, (2) redefi- 
nition of the centerline space curve to satisfy design constraints, and (3) perturbations of the 
cross-sectional shape to reflect specified design iterations. The approach taken here is to develop 
a geometry' pre-processor to augment the existing geometry and grid generation programs for 
internal inlet duct configurations. 

In solving the equations of fluid dynamics by numerical means, two major problem areas 
are encountered: the generation of the three dimensional geometry and grid, and the algorithm 
by which the governing equations are processed to obtain a solution. Grid generation envolves 
three processes: (1) description of the boundary surface, (2) construction of the grid lattice within 
the boundary surface, and (3) solution dependent adaptation or clustering of the grid lattice to 
accurately resolve regions of rapid solution variations, i.e. large shear. 

Within this framework, the modern computational codes have separated the functions of 
grid generation and flow solution, i.e. the geometry and mesh construction has been divorced 
from the solver. This separating of the functions is a sensible choice within a fundamental re- 
search environment where the processes of grid generation and flow solution are destinctly dif- 
ferent and where the design is fixed, i.e., where variations of the same geometric configuration 
need not be considered. However, in a analysis-design engineering environment, or even in 
analysis of many three dimensional viscous flows, large advantages can be gained by having the 
ability to modify an existing gridfile without returning to a separate geometry and mesh genera- 
tor. For example, adaptive mesh clustering is a very' powerful technique to effectively use a lim- 
ited number of grid points to define the multiple length scales that can arise in fluid dynamics 
problems, such as boundary layer and inviscid core regions. In numerous problems, particularly 
in three dimensional space, the location of the regions requiring definition are not known a-priori, 
but emerge as the solution develops. A vortex ingestion-boundary layer interaction is a case in 
point. For such problems, adaptive mesh re-clustering of an existing gridfile, without the neces- 
sity of returning to a separate mesh generator, becomes a very powerful computational technique. 

Redefinition of the centerline space curve artd perturbation of the cross-sectional shape of 
the existing gridfile inlet duct description can have important significance within a analysis-design 
environment, where analysis of many variations of the same basic geometry with different inflow 
conditions is often the design task. Special constraints can easily be imposed on the inlet duct 
geometry within the system of patitioned geometry and mesh generation. For example, the 
centerline definition can easily be altered without affecting the inlet flow area distribution. This 
is shown in Figure (2), where an inlet duct equivalent to the original F/A-18 configuration was 
constructed with a straight rather than an out-of-plane offset centerline. Likewise, the cross- 
sectional shape can easily be modified, without affecting either the centerline space curve or flow 
area distibution. This also is shown in figure (2), where an equivalent circular cross-section 
F/A-18 inlet duct is presented. These are very powerful computational capabilities within an 
design environment. 


Reference Coordinate System 

In the present development, the primary flow direction and body fitted coordinated system 
are related to a smooth reference line associated with the geometry of interest. This reference 
line, usually chosen to be the duct centerline, is a space curve required to have continuous second 
derivatives. The planar surfaces normal to the reference line at each axial location intersect the 
duct surface and are identified as the transverse (cross-sectional) planes of the duct geometry. 
The primary flow direction is taken as perpendicular to the transverse plane at each axial lo- 
cation, and thus parallel to the reference line. The reference line thus provides a necessary link 
between the duct geometry and the flow approximations made to permit solution as an inial value 
problem. An orthogonal reference coordinate system is derived to fit the reference line. In this 
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system, the transverse coordinate are "Cartesian" and lie in the transverse plane, and the axial 
coordinate is normal to the transverse plane. Although the orthogonal coordinate system does 
not fit the duct geometry, it is aligned with the primary flow direction and thus provides a con- 
venient method for introducing the flow approximations in a simple form. The approximate 
equations are written in this orothogonal reference coordinate system. A nonorthogonal body 
fitted coordinate is then devised by transforming the transverse Cartesian coordinates to obtain 
new transverse coordinates which fit the intersection of the duct with the transverse normal 
planes, and by simultaneously transforming the axial coordinate to fit the duct geometry. Tran- 
sformations from the body-fitted coordinate, to the computational orthogonal coordinates, x„ 
built on the reference line are calculated as part of the body-fitted system. The governing 
equations are then solved following a Jacobian transformation to the body fitted coordinate sys- 
d yj 

tern, i.e. Thus, for example: 


dx, 


du n 


dyj 


du„ 


dxj dx t dy/j 


( 1 . 1 ) 


and in general 


d _ ty d 

dx i dx t dyj 


( 1 . 2 ) 


The Governing Equations 

The reduced form of the Navier Stokes equations are derived through approximations 
made relative to the orothogonal coordinate system built upon the reference line described in the 
previous section. The reference line specification links the flow approximations with the duct 
geometry. Transverse coordinate surfaces are planes normal to the reference line. Equations 
governing the primary flow velocity, U p , and a secondary vorticity, co normal to the transverse 
coordinate surfaces are derived utilizing approximations which permit solution of the equations 
as an initial value problem. Terms representing diffusion normal to the transverse coordinate 
planes are neglected. Approximate pressure gradients are imposed in the streamwise momentum 
equation from an apriori known pressure field obtained from a solution of the potential flow in 
the duct. Secondary flow velocities are determined from scalar and vector potential calculations 
in transverse coordinate surfaces, once the primary velocity is known. 

The analysis is based on the decomposition of the overall velocity vector field U into a 
primary flow velocity U p and a secondary flow velocity U s . The overall or composite velocity is 
determined from the superposition 

U = U p +U s (2.1) 

The primary flow velocity is represented as 

U p = u p i n (2.2) 

where i n is the normal to the transverse plane. The primary flow velocity u p is determined from 
the solution of the primary momentum equation. In general orthogonal coordinates, the primary 
momentum equation becomes 
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pu p du p 
h x dx, 
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pv s du p 
h 2 dx 2 
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+ ( JRe ) 

+ (JRe) 
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pw s du P 

h 3 dx 3 


d ^1 ^3 d 
P 


+ u P ( v s^n + W X\ 13 ) ~ V *K 21 — w j^ 3 i ~ 


h x dx x 


(Pj+P v ) 


dx 2 h 2 


d M2 
dx 3 h 3 


dx^ 


dx* 


jr) + ^ K --t 


+ *?***> -t 


(2.3) 


dh, 


where J = h,h 2 h 3 and K 0 = (h,hj)~ l 

The streamwise pressure gradients are approximated from an inviscid flow obtained by an 
a priori analysis, and are corrected for viscous effects by a mean pressure gradient term which 
depends only on the streamwise coordinate. In equation (2.3), P, is an a priori known pressure 
field, and P v is the one-dimensional correction to account for the mean viscous pressure loss. 
For internal flows, P v is determined to ensure that the integral mass flux condition is satisfied: 


m 


= f pUpdA = CONSTANT 

J A 


(2.4) 


The secondary velocity U s is derived from scalar and vector surface potentials denoted as 
(j) and i j/, respectively. For p the density and p 0 a constant reference density, U s is written as 


where V, is the surface gradient operator defined by 

v, = v-/ n (?,.v)' 

The scalar and vector potential equations become 


d 

M3 

d<f> + 

d 

M 2 

d4> 

d 

dx 2 

P h 2 

dx 2 

dx 3 

P h 3 

dx 3 ~ 

dx- 

9 ( 

h 1 

dh x \j/ 

V 

9 (- 

h 2 ! 

dh x >jj 

dx 2 { 

h x h 2 P 

dx 2 

r 

dx 3 \ 

h x h 3 P 

dx 3 


— h 2 h 3 ^2 


and the vorticity transport equation in general orthogonal coordinates can be written 


(2.5) 

( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


P P 9 ■ h 2 h 3 a + h 3 pv£l + -53- h 2 pw^l 


dx. 


dx - 1 


+ 

+ 


Sh 3 w s g 
dx x dx 2 

dpVs dv s 

dx 2 dx 3 


P“» 


dh-,v, 


■2 v s 


P“p 


dx, dx 1 


dx-, 


(h 3 pUpK x3 ) + {h 2 pu p K x2 ) 


dpv s dv s dpw s dw s dpw s dw s 


dx 3 dx 2 


= Re~ l 


+ ■ 


dx 2 dx 3 
d h 2 d 


dx 2 h x h 2 dx 2 ' dx 3 h x h 3 dx 3 


dx 3 dx 2 


(ph x Q) 


(2.9) 


where the decomposed transverse velocity components are given in orthogonal coordinates by 
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( 2 . 10 ) 


v * “ v + v 


1 + 1 


1 dh^ 


W^W'fi + W'j, 


h x dx 2 h x h 3 P dx 3 
1 d<j> ii dh x \l> 


h 3 dx 3 


h x h 2 P dx 2 


The equation of state for a perfect gas is given by 

p = pRT 


(2.11) 


(2.12) 


where R is the gas constant and T is temperature. In the present study, the stagnation energy 
was assumed to be a known constant, E 0 , and the energy equation can be omitted from consid- 
eration. For constant stagnation energy, the gas law can be written as 


Pi+Pv~ 

where y is the specific heat ratio. 


1 ( r U'U \ 

2 ) 


(2.13) 


The Turbulence Model 


The mixing length model used in this analysis employes an an eddy-viscosity formulation 
for the Reynolds stresses expressed in the form 



P-t d'J 
Re dxj 


(3.1) 


where the effective viscosity p T is added to the laminar viscosity p. The turbulent viscosity is re- 
lated to the mean flow variables by means of a mixing length distribution 


^ = p/ 2 (2^:I) 1/2 (3.2) 

where e is the mean flow rate of strain tensor defined by 

^ = y[(Vv) + (ViT) 7 ] (3.3) 

The mixing length € is determined from the empherical relationship of McDonald and Camarata, 
Ref. 11, for equilibrium turbulent boundary layers, which can be written as 


€(y) = 0.09<5 6 tan h 




0.093, 


D 


(3.4) 


where 5 k is the local boundary layer thickness, k is the von Karman constant, andj' is the distance 
from the wall, and D is the sublayer damping factor defined by 

n-r^zll 

D v p <7 j (3-5) 


where 


y 


+ 



(3.6) 
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and p is the normal probability, t is the local shear stress, y* = 23, ct, = 8.0 , andy is the distance 
to the wall. 


The FLARE Approximation 

The analysis as presented here is applicable only when the primary velocity is not negative. 
Since "small" regions of reverse flow can arise in curved inlet ducts, the numerical method is lo- 
cally modified to permit forward marching when the flow contains small regions of reverse flow. 
The technique used follows Reyhner and Flugge-Lotz, (Ref. 12), by adding small artificial con- 
vection at grid points where the primary flow is reversed. This is known as the FLARE approx- 
imation, after the authors of Reference 12. For thin regions of reverse flow, although the area 
of flow separation can be very large compared to the passage itself, the technique permits the 
analysis to proceed downstream beyond reattachment, confining the FLARE approximation to 
the separated region. 


Vortex Generator Model 


One of the most commonly used methods to control local boundary layer separation en- 
tails the placement of vortex generators upstream of the problem area. Vortex generators in use 
today are small wing sections mounted on the inside surface of the inlet duct or wing surface, 
inclined at an angle to the oncoming flow to generate a shed vortex. The generators are usually 
sized to local boundary layer height to allow for the best interaction between the shed vortex and 
boundary layer itself, and are usually placed in groups of two or more upstream of the problem 
area. The principle of boundary layer control by vortex generators relies on induced mixing be- 
tween the external or core stream and the boundary layer region. This mixing is promoted by 
vorticies trailing longitudinally over the duct surface adjacent to the edge of the boundary layer. 
Fluid particles with high momentum in the stream direction are swepted along helical paths to- 
wards the duct surface to mix with and, to some extent, replace the low momentum boundary 
layer flow. This is a continuous process that provides a source of re-energization to counter the 
natural boundary layer growth caused by friction, adverse pressure gradients, and low energy 
secondary flow accumulation. There are two basic configurations of vortex generators. In one 
configuration, all the vortex generators are inclined at the same angle with respect to the 
oncoming flow direction. These are called co-rotating configurations because the shed vorticies 
rotate in the same direction. In the other configuration, the vortex generators are grouped in 
pairs inclined in the opposite direction to the flow, such that pairs of counter rotating shed 
vorticies are generated. 

Co-rotating vortex generators are very competitive in reducing flow separation if the gen- 
erators are properily selected and located. The main advantage of co-rotating type vortex gen- 
erators are their downstream effectiveness resulting in more effective usage of the vortex energy 
within the affected boundary layer. According to design "wisdom", this type of vortex has a few 
special advantages when used within S-duct inlet configurations, namely: (1) the induced vortices 
will remain close to the wall; consequently a "cleaner" core flow will result, and (2) the induced 
vorticies will counteract the natural and often strong secondary flows which develop. 

Counter-rotating, equal strength vortex generators have been used in a number of aircraft 
inlet ducts, the F/A-18 and the center inlet duct on the production 727 aircrafts, to name a few 
examples. This type of vortex generator is very effective in reducing flow separation if the vortex 
generators are placed slightly upstream of the region of separation. According to design "wis- 
dom", the disadvantages of this type of generators, as compared to co-rotating generators, are 
that the induced vorticies tend to lift off the duct surface, thus reducing their effectiveness, caus- 
ing higher loss in inlet recover)' and larger total pressure distortion at the engine face. 

The performance of vane-type vortex generators was evaluated by Taylor, (Ref. 13), for 
diffusers and airfoils at low speed, and by Valentine and Carrol, (Refs. 14 and 15), for airfoils and 
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wings at high speeds. This work provides trends in effectiveness for certain vortex generator de- 
sign variables such as angle-of-attack, height, distance ahead of separation, etc. Attention was, 
however, focused on the detailed changes that were produced in the boundary layer as a result 
of placement of vortex generators in the flow. Percy and Stuart, (Ref. 16), extended the study 
of the effects of various design parameters and concluded that the strength and disposition of the 
individual induced vortices was more important than the details of the boundary layer upstream 
of the imposed pressure gradient. These results suggest that modeling of the vortex generators 
within a marching RNS analysis could be be very successful, probably because of the relative 
length scales encountered between the vortex generator region and the inlet duct itself. 

The model for the vortex generators within the RNS analysis, described by Kunik, (Ref. 
17), takes advantage of the stream function-vorticity formulation of the governing equations. 
The shed vortex is modeled by introducing a source term into the vorticity equation that is a 
function of the geometric characteristics of the generators themselves. This source term is in- 
troduced at every point in the cross-plane in the form of the following expression 

r p = i>-( Ci ' 2 ) (4.1) 

where T, is the vortex strength at any point in the crossplane, T 0 is the vortex strength at the tip 
of the generator, r is the distance between the field point and at the tip of the generator, and c, 
is a constant which controls the decay of the shed vortex strength in the crossplane. This vortex 
model resembles the one proposed by Squire, Ref. 18, except that it neglects the variation of 
viscosity in the crossplane. 


Inlet Duct Geometry and Flow Descriptors 


The reference line or centerline, in addition to defining a coordinate system in which to 
make the calculations, also characterizes the inlet duct in terms of a parameter t. Thus, the duct 
contours can be represented in terms of the centerline space curve 

X c1 = X{t) (5.1) 

Y cl =Y( t) (5.2) 

Z c i = Z(t) (5.3) 

with prescribed cross-sectional planes normal to the centerline with flow area A( t), and a ge- 
ometric shape factor SF(t) defined by 


SF = 1.0 



(5.4) 


which characterizes the cross-sectional shape of the duct. In equation (5.4) D hyd is the local hy- 
draulic diameter, A the local flow area, and r{6) is the local cross- sectional radius of the inlet duct 
in the local transverse plane. As the cross sectional shape approaches a circle, SF will approach 
the value of 1 .0. 


There are a number of flow discriptors that will be use to characterized the effects of inflow 
distortion and its transmission through the inlet duct. These include the static pressure coefficient 
and the mass flow averaged static pressure coefficient defined by 


Cp = 


(P-Po) 

% 


(5.5) 
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(5.6) 



where m is the mass flow given by 


m -= 


l 


pu p dA 


(5.7) 


and p 0 and q 0 are the reference static and dynamic pressure. Likewise, the total pressure coeffi- 
cient and the mass flow average total pressure are defined by 


CP: 


(Pt~P o) 

% 


(5.8) 


Pt=±^J>U p {Pt)dA (5.9) 

Thus, the mass flow averaged inlet total pressure recovery n d is given by 

P^ef 

* d = -=- (5. 1 0) 

where the subscript ( ef) refer to the engine face station or aerodynamic interface plane, and (/) 
indicates the inflow station. 


In order to quantify effects of vortex disturbances entering the inlet, a number of param- 
eters will be defined to characterized this type of disturbance. These parameters include the sec- 
ondary flow kinetic energy 


KE S 




(5.11) 


the local swirl kinetic energy 



and the local primary flow kinetic energy 



(5.12) 


(5.13) 


The ratio of local secondary to primary flow kinetic energy becomes 

KE S 

k *=~ke; 

and the ratio of swirl to primary flow kinetic energy is given by 

ke q 

K ° KE p 


(5.14) 


(5.15) 
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Steady State Engine Face Distortion Descriptors 


It is impractical to measure anything at the engine face when the engine is installed and 
operating, consequently, the engine and inlet designers agreed upon an Aerodynamic Interface 
Plane (AIP) which is forward of the compressor face but sufficiently close to the engine face to 
have a similar flow field. Current U.S. practice uses forty or forty eight transducer probes ar- 
ranged in eight rakes with five or six rings. The radius of each ring is set such that all probes are 
at the centroid of equal areas. All distorton descriptors, whether they quantify steady state or 
transient distortion conditions, are always calculated relative to the standard rake located at the 
AIP. 


The most widespread quantitative distortion descriptor available in the literature, because 
of its use in the earliest measurements on inlet ducts in the late 1950's, is simply: 

Dt = ( ~ Plmin ) (6.1) 

\ **ave } 

where Pt max is the maximum rake total pressure, Pt mm is the minimum rake total pressure, and 
Pt m , is the area weighted average rake total pressure. In experimental data reduction, it is as- 
sumed that the both the static pressure and temperature are constant and steady across the 
aerodynamic interface plane (AIP); thus both the velocity and Mach number can be considered 
functions only of total pressue and the distribution of this quantity is the only measurement that 
needs to be made. This parameter is always useful to determine for comparison purposes and to 
describe the 'general health' of inlet ducts irrespective of the type of powerplant that may be used. 

More advanced distortion descriptors, introduced in the late 1960's and 1970's, take into 
account the Dt distortion of each ring of total pressure measurements. Thus, the radial distortion 
Dt rad is defined as 

rj f ^max ~ P^ave \ ff - ~ 

Dt rad—\ pi /ring ( 6 - 2 ) 

where Pt m , is the average total pressure for a given ring radius and Pt m „ is the maximun local ring 
total pressure. The circumferential distortion Dt cir is defined as: 

(63) 

where Pt min is the lowest average total pressure in any 60-degree or 180-degree arc for a given ring 
radius having an average ring total pressure Pt m ,. 
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RESULTS AND DISCUSSIONS 
Comparison and Validation with Experimental Data 

To demonstrate the accuracy of the numerical procedure obtained with the 3D RNS code 
for internal duct flows typical of high angle-of-attack conditions, a series of numerical simulations 
were carried out on the University of Tennessee diffusing S-duct validation case. Vakili, Wu, 
Liver, and Bhat, Ref. 10, in this experimental investigation sponsored by NASA Lewis Research 
Center, obtained a series of measurements in a 30-30 degree diffusing S-duct of area ratio 1.5 with 
and without vortex generators, Figure 3. The 30-30 degree circular cross-section S-duct, shown 
in Figure 3, was made from two symmetric sections. The inside hydraulic duct diameter D, was 
16.51 cm. and the mean radius of curvature R was 83.82 cm. A 72.20 cm. straight pipe was in- 
stalled upstream of the curved section to allow for the development of the turbulent boundary 
layer to the desired thickness. Another pipe of length 152.40 cm. was installed downstream of the 
S-duct. 

All measurements were made at an inlet Mach number of 0.60 at the reference measure- 
ment station in the straight section X= 24.8 cm. The flow parameters at this station were used 
as reference conditions for non-dimensionalizing the experimental data. The experimental survey 
stations correspond to X/Z), = 0.0, 1.29, 2.49, and 5.2. At each survey station, a five port cone 
probe was traversed radially at ten azimuthal angles, approximately 20 degrees appart, on both 
sides of the symmetry plane. At least seventy points were measured at each traverse. 

A polar grid topology, (Figure 4), was chosen for the analysis of the University of 
Tennessee diffusing S-duct, consisting of 49 radial, 49 circumferential, and 101 streamwise nodal 
points in the half-plane. The internal grid was constructed such that the transverse computa- 
tional plane was perpendicular to the duct centerline. Grid clustering was used both in the radial 
and circumferential directions by redistributing the nodal points along these coordinate lines to 
resolve the high shear region near the wall and the separation region in the second bend. The 
flow in the inlet duct was turbulent, with an entrance Mach number of 0.6, Reynolds number 
based on hydraulic diameter of 1.76.rl0 6 , and the inflow corresponds to a shear layer approxi- 
mately <5/Z), = 0.05. These initial conditions were applied at approximately an axial station 1.54 
inlet diameters (D,) upstream of the duct entrance. 

Figure 5 shows a comparison between the experimental and computed total pressure co- 
efficient contours at XjD, = 5.2 . In both the experiment and analysis, the flow in the S-duct 
separated and reattached in the second bend upstream of XjD t = 5.2. This was due both to the 
adverse pressure gradients and the effect of pressure-driven secondary flow resulting from duct 
curvature. Experimental measurements and computational results from the 3D RNS code in- 
clude the phenomena of separation and reattachment upstream of X/Z), = 5.2 and show excellent 
agreement for a simple mixing length turbulence model. 

As separation was encountered in the second bend of the S-duct, three pairs of vortex 
generator devices were installed in the duct at X/D t = 0.09, and at azimuthal angles of -38.0, 0.0, 
and 38.0 degrees relative to the top of the inlet, and with respect to the streamwise direction. 
The vortex generator pairs had geometric incidence angles of +16.0° and —16.0° relative to the 
duct centerline. Figure 6 shows the comparison between the experimental and computed total 
pressure coefficient contours at X/Z), = 5.2. Comparison of contour levels between the separated 
case, Figure 5, and the vortex generator case, Figure 6, show that the vortex generators success- 
fully mixed the high energy core flow with the low energy flow in the wall region to supress sep- 
aration. In general, the interaction between the induced vortex generator flow and the pressure 
driven secondary flow that was computed was physically realistic and the agreement between 
experiment and analysis is considered very good. 
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Figures 7 and 8 show additional flow characteristics obtained with the 3D RNS analysis 
with the vortex generator modeling. The secondary flow structure from the vortex generator 
model just downstream of the generator region, i.e. at XID, = 0.18 is shown in Figure 7, and 
clearly reveals the three pairs ofvorticies that arise from the three pairs of counter-rotating gen- 
erators. The limiting streamline signature shown in Figure 8 indicates that the generator config- 
uration tested elliminated the flow separation encountered in the second bend, and reveals the 
familiar topographical pattern through the generator region itself. 

Flow Separation and Vortex Liftoff 

The three dimensional separation encountered in the University of Tennesse diffusing S- 
duct was very large in area, Figures 9 and 10, but thin in the direction normal to the wall such 
that it did not alter the pressure distribution in a substantial manner. Secondary flow resulting 
from duct curvature caused an accumulation of boundary layer near the innerwall of the first 30 
degree bend. The thick boundary layer thus established was especially susceptible to flow sepa- 
ration because of the adverse pressure gradients induced by the reverse curvature section of the 
second 30 degree bend. A comparison between the computed oil flow patterns (represented by 
the limiting streamline topology in figure 9) and the experimental oil flow patterns presented in 
figure 10 shows excellent correspondence. Of importance is the fact that the viscous marching 
RNS analysis method using FLARE approximations captures the reverse flow region of this 
separation in a single pass calculation. 

To meet the required confidence level of code validation, it is also essential that the 3D 
RNS marching analysis be able to capture the known topological structure of the limiting 
streamlines in the vicinity of separation. The very important and striking features of the sym- 
metric pair of spiral nodes and saddle points were clearly captured by the Reduced Navier Stokes 
analysis. This very familiar topological pattern is known to describe the important stage in the 
development of the pair of counter rotating vorticies that form in the first 30 degrees of turning 
that results in vortex liftoff in the second bend, Figures 9 and 10. Vortex liftoff and the resulting 
separation are jointly a major cause of total pressure loss and distortion at the engine face. 

In general, the flow separation studied here was dominated by pressure forces rather than 
shear forces, as would be the case with massive separation in the inlet duct, such that the influ- 
ence of the stress-driven flows was small and the effective viscosity approach surprisingly suc- 
cessful. These separation interactions could be called "vorticity separation" since they are 
dominated primarily by the transport of vorticity through the inlet duct (see following section). 
The case with vortex generators was also dominated by pressure driven secondary flows, as tes- 
tified by the very good results from a simple modeling of the vortex generators in terms of a 
vorticity signature calculated from the physical properties of the generators themselves. The class 
of problems dominated by pressure forces rather than shear forces are known as vorticity domi- 
nated flows, and make up a very important class of internal flow problems. 

Vorticity Dominated Internal Flow Fields 


A predominant direction of flow can be identified in many high Reynolds number flow 
problems, both extemaly and intemaly. This primary direction may be derived from the direction 
of flight in external aerodynamics or from the ducted nature of internal flows. Viscous effects 
are very important, but they are generally limited to the region near the wall. In curved ducts, 
the outer portions of the boundary layer away from the duct walls behave as a rotational inviscid 
flow region in which secondary flow is generated (in the terminology of secondary flow theory) 
by turning of the transverse shear. This rotational inviscid behavior is reminiscent of the middle 
rotational inviscid layer of triple-deck theory for a two-dimensional boundary layer subjected to 
a streamwise disturbance, although the latter case does not involve streamwise vorticity. As the 
flow proceeds downstream, the secondary flow convects the streamwise vorticity and distorts the 
primary flow in what amounts to a strong interaction in the transverse directions between viscous 
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and inviscid (rotational or potential) flow regions. Significant regions of rotational but essentially 
inviscid flow occur as a result of this process, and the potential-core region will eventually dis- 
appear. 

Several numerical and physical factors can affect accurate predictions of high Reynolds 
number vorticity-dominated internal flows. These include the effects of numerical smoothing, (in 
the case of Full Navier Stokes analyses), the turbulence model, and the need for sufficient fine 
grids to resolve the details of the viscous boundary layer in the near wall region where the 
vorticity is the largest. For the computations of the flow within the University of Tennessee dif- 
fusing S-duct, two calculations were made to study the effects of near wall grid resolution. In the 
first case, the radial grid spacing was chosen to give a nominal y+ value of 8.5 at the first grid 
point above the duct wall, and in the second case, a nominal y + value of 0.5 was chosen. Figures 
1 1 and 12 show the effects of this near wall grid lattice resolution on the total pressure coefficient 
contours, (Figure 11), and the secondary flow structure, (Figure 12), both at the engine face sta- 
tion, X/Z), = 5.2. The very strong effect of near wall grid resolution on the structure and strength 
of the engine face flow results from the very nature of vorticity dominated internal flows. Since 
the largest value of vorticity occurs in the near wall region, and secondary flow is generated by 
turning of streamwise vorticity (shear), a very strong effect of near wall grid resolution was thus 
realized as stronger secondary flow (Figure 12), which had consequently an appreciable influence 
on the primary flow (Figure 11). This is really an inviscid rotational phenomenon rather than 
a viscous interaction, and indicates that secondary flow has its ultimate origin very near the wall. 
The influence of the turbulence model will thus also strongly influence the strength of secondary 
flow, but only as an invsicid rotational phenomena through the near wall vorticity distribution 
and not necessarily as part of the turbulent properties of the flow. 

Figure 13 presents the limiting streamline signature associated with vortex liftoff obtained 
from the RNS code for the y+ = 0.5 calculation, while Figure 14 shows the partical traces of 
vortex liftoff, again for the y + = 0.5, calculation for particles originating at>>/<5, = 0-03 It is quite 
apparent that fluid particles that originate well inside the entrance boundary layer can influence 
the core region of the engine face station. 

Fundamental Internal Vortex Interactions 


Unlike external aerodynamics where the two-dimensional flow assumption is valid over 
extended regions of the flow field, the internal flow problem is plagued with three-dimensional 
flow effects. When viscous effects or other sources of vorticity are present, three-dimensional 
flows differ fundamentally from their two-dimensional counterparts in that large secondary flows 
are generated by a deflection of the primary flow and/or other mechanisms. Secondary flow 
theory (reviewed by Horlock and Lakshminaryana, (Ref. 20), and Lakshminarayana and Horlock, 
(Ref. 21) affords considerable insight into the generation of secondary flow and establishes that 
large secondary flows can be generated by small deflection of vorticity or shear. The large sec- 
ondary flows thus generated often exert an appreciable influence on the primary flow, and thus 
aerodynamic performance, viscous losses, and engine face distortion can be significantly affected. 

The geometry of the bounding walls and flow area distribution are the primary drivers for 
3D flow effects in internal passages. The development of secondary flows in transition ducts or 
curved ducts provide examples of 3D flows generated internally. The effect of the three- 
dimensional flow field on vorticity dynamics is described by the phenomenon of vortex stretching 
(reviewed by Batchelor, Ref. 22). The term co • VU in the _vorticity dynamics equation, (see for 
example, reference 22), where <y is the vorticity vector and U is the velocity vector, describes the 
time rate of change of vorticity due to "stretching" of the vortex filaments by the 3D flow field. 
This phenomenon, i.e. vortex stretching, is purely inviscid in origin. 

Another fundamental internal interaction that can arise is when vorticity is introduced 
externally, i.e. either through a vortex ingestion or generated internally by means of vortex gen- 
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erators. The resulting flow that develops under these circumstances has its origin in inviscid flow 
theory as well. Batchelor, (Ref. 22), describes the spiral motion of an initially straight streamwise 
vortex filament in a cylindrical duct with area contraction. The qualitative nature of this flow 
was such that as the vortex passed through the contraction section, it moved radially inward, and 
the azimuthal velocity of a material point on the vortex-line increased according to the rule 
rco = const, hence it increased. Therefore, the predominant characteristic of vortex motion within 
an inclosed duct is the spiraling of the vortex core in a direction having the same sense as the 
vortex. Viscous effects are important, but their influence is felt by increasing boundary layer 
blockage thus causing the vortex-filament to move radial inward. The azimuthal velocity will 
increase as a result of rco = const, behavior but decrease as a result of diffusion of vorticity from 
the vortex axis. 


Vortex Interactions within the F/A-18 Inlet Duct 

Two types of vortex-boundary layer interactions will be examined in this study. In the first 
case, vorticity is generated internally by means of one pair of counter-rotating vortex generators 
within the F/A-18 inlet duct. The purpose of these generators was to suppress flow separation 
that was likely to occure downstream of the a small "bump" within the duct, rather than to con- 
trol engine face distortion. The second type of vortex-boundary layer interaction that was ex- 
amined was under the condition where vorticity was externally introduced by means of an 
ingested vortex. 

The F/A-18 inlet duct geometry definition used in this study is is shown in Figure 15, and 
the single block polar grid topology, with 49 radial, 98 circumferential, and 101 streamwise points 
within this inlet duct is presented in Figure 16. The computations were made at an inlet entrance 
Mach number of 0.6, and Reynolds number of 8 -OjcIO 5 based on the hydraulic inlet diameter (Z),) 
, and inflow conditions that correspond to a shear layer thickness (5/Z), = 0.056. 

The vortex generators used in this study were located at a streamwise station X/Z) f = 3.53 
upstream from the engine, and at an azmuthhal angle of 180.0 degrees relative to the inlet vertical 
axis, (Figure 17). All vortex generator dimensions, (Figure 17) are normalized in terms of the 
hydraulic inlet diameter (Z)<), which was 23.34 inches. The vortex generator incidence angles were 
+ 16.0 and -16.0 degrees with respect to the axial direction. Shown in Figures 18 and 19 are the 
limiting streamline signatures without and with vortex generators respectively. The characteristic 
vortex signature of the limiting streamline is through and downstream of the vortex generator 
reqion quite evident in Figure 19. An interesting study that should be made is to determine the 
effect of vortex generator flow alignment on aerodynamic effectiveness. 

The second type of vortex-boundary layer interaction to be examined is under the condi- 
tion where vorticity is introduced by means of a vortex ingestion. The nomenclature used for the 
constuction of the ingested vortex is shown in Figure 20. Generic entrance conditions were es- 
tablished in each case by centering the core of the vortex of dimension (r co JR) at (rJR„ 6 C ), where 
(/?,) is the hydraulic inlet radius. It must be understood that the detailed inflow conditions of a 
real vortex ingestion is largly unknown and must be determined experimentally, however, this 
study seeks only to represent the starting inflow conditions in a generic sense such that studies 
can be made to relate inlet performance, i.e. inlet recovery, and engine face distortion, to the ge- 
ometric properties of the duct. Tables I and II present the range of test cases considered for this 
study, which have been arranged in three groups. In the first group, (Table I), the vortex entry 
point was held fixed at an (rJR it 9 C ) of 0.779 and 225.0° while the swirl angle decreased from 0.0 
to —45.0°, i.e. (counter-clockwise). The second group of cases, (Table II), also kept the entry 
point fixed at an ( rJR it 6 C ) of 0.779 and 225.0° , while the swirl angle increased from 0.0 to 
40.0°, i.e. (clockwise). In the third group of cases, the swirl angle was held fixed at —30.0° , while 
the radial location of the vortex core (r e /Z?,) was increased from 0.657 to 0.987. 
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Figures 21 and 22 show the inflow vortex structure at the entrance, (Figure 21), of the inlet 
duct X/A=0.0, and that vortex structure which would be experienced at the engine face station, 
(Figure 22), at X/Z),= 6.5, for both the clockwise vortex and counter-clockwise vortex, with sw'irl 
angles of +10.0° and -10.0° , respectively. The RNS solution within the F/A-18 inlet duct 
showing the limiting streamlines with an ingestion of a +10.0° clockwise vortex is shown in Fig- 
ure 23, while the limiting streamline pattern produced by the of -10.0° counter-clockwise vortex 
is shown in Figure 24. In general, the entrance vortex wants to spiralal around the inside surface 
of the inlet duct in a direction corresponding to the sense of the initial vortex, i.e. whether the 
vortex is rotating clockwise or counter clockwise. This spiralling characteristic around the inside 
wall of the inlet duct is indicated in Figures 23 and 24 by the coalescence of the limiting stream- 
lines. This coalescences of the streamline reveals the vortex surface trajectory, and indicates a 
secondary (transverse) flow separation within the duct. The influence of the vortex generator can 
be seen for the clockwise vortex ingestion, (Figure 24), but interestingly, the counter clockwise 
vortex impinges on the vortex generators, deflecting the vortex trajectory, and rendering the 
generators ineffective. 

In order to quantify the effects of a vortex disturbance entering the inlet duct, both the 
secondary to primary flow ( k s ) and swirl to primary flow (k # ) kinetic energy ratio were defined 
by equations (5.14) and (5.15) respectively. Figures (25) and (26) presents both the swirl and 
secondary flow the kinetic energy signature of the +10.0° clockwise vortex ingestion, (Figure 25), 
and —10.0° counter-clockwise vortex ingestion, (Figure 26), both compared to the background 
signature, i.e. a, = 0.0°. It is assumed that the swirl kinetic energy component (k„) is that portion 
of the secondary flow kinetic ratio (k s ) that would be elliminated with inlet guide vanes, (IGV's). 
IGV's are present in all American modern military aircraft high performance engines except the 
TF41, which is a derivative of the European Spey. 

Thus, these parameters provide a unique signature which can track and quantify the pas- 
sage of each and every vortex disturbance entering the inlet duct. Therefore, a relationship be- 
tween the inlet duct geometry, described in terms of a distribution of area, A( r), and shape factor, 
SF(t), along a centerline space curve with geometric properties of slope and curvature, can be 
established. This provides a quantitative methodology whereby the aerodynamic flow quantities 
characterizing the vortex ingestion can be directly related to the important geometric quantities 
characterizing the inlet duct geometry. 

Figures 27 through 30 present the particle traces released from y/t 5, = 1.14 for vortex 
ingestion swirl angles a, = 0.0°, -10.0° , -20.0°, and -30.0°. The spiralling characteristics of the 
ingested vortex are clearly revealed in these figures along with the characteristic that as the initial 
swirl angle increases, the vortex helical angle increases, i.e. the vortex spiralis further around the 
inside surface of the duct. Particle traces resulting from an ingestion of a counter-clockwise 
vortex of strength a, = —40.0°, (Figure 31), reveal very clearly the helical-like vortex trajectory 
as seen from a downstream view. Figure 32 presents the trajectory of the vortex center relative 
to the inlet duct centerline for the a, — -40.0° vortex ingestion case. 

The effect of vortex strength on the engine face impingement location is shown in Figure 
33. The a, = 0.0° impingement location represents the ingested vortex entry point at the inlet 
face. Thus, as the strength of the ingested vortex increases, its compressor face angular 
impingement location increases and the core of the vortex moves radially inward. The radial in- 
ward movement is a results of increased boundary layer blockage caused by the increased strength 
of the ingested vortex. Figure 34 presents the effect of radial entry point on engine face vortex 
impingement location, i.e. the third set of vortex ingestion test cases listed in Table II. At a 
vortex entry position greater than r c /A = 0.779, the inflow vortex is "cut-ofF by the lip, and 
consequently, its strength decreases, i.e. the circulation around any closed path surrounding the 
vortex filament decreases. As a consequence, a maximum angular impingement location occurs 
at the engine face station, (Figure 34). 
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The overall qualitative nature of a vortex ingestion, at the Reynolds number considered, 
was basically "inviscid like", i.e. the dominant aerodynamic characteristics have their origin in 
inviscid flow theory. Therefore, the predominant characteristic of vortex motion within an in- 
closed duct is the spiralling of the vortex core in a direction having the same sense as the vortex. 
Viscous effects are important, but its influence is felt by increasing boundary layer blockage, thus 
causing the vortex-filament to move radial inward. The azimuthal velocity will increase as a re- 
sult of rto = const, behavior but decrease as a result of diffusion of vorticity from the vortex axis. 
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CONCLUSIONS 


The goal of the present investigation was to combine the Reduced Navier Stokes solution 
technique with the concept of partitioned geometry and mesh generation to form an efficient 3D 
RNS code aimed at the CFD analysis-design engineering enviornment. Partitioned geometry and 
mesh generation is a geometry pre-processor to augment existing geometry and grid generation 
programs and has the advantage that the solver can (1) recluster an existing gridfile mesh lattice, 
and (2) perturb an existing gridfile geometry and mesh to alter the cross-sectional shape or 
centerline definition without returning to the external geometry and grid generator. 

The present results provide a quantitative validation of the initial value space-marching 
3D RNS procedure and demonstrates accurate predictions of the engine face flow field, with a 
separation present in the inlet duct as well as when vortex generators are installed to supress 
separation. The computing time, CPU = 6.5 min. on the CRAY XMP for 2.28 x 10 s grid points, 
for both the baseline case and the case with three pairs of counter-rotating vortex generators, is 
sufficiently rapid for routine use in an analysis-design engineering environment. 

Initial value space marching 3D RNS procedures using FLARE approximations can 
adequatly describe the topographical features of 3D flow separation associated with vortex liftoff 
within inlet ducts. The success of this RNS analysis in describing this phenomenon is due to the 
existence of a class of separated flows, which can be called "vorticity separations'', which are 
dominated by the transport of vorticity rather than turbulent shear effects. 

Adequate resolution of turbulent flows including viscous sublayer resolution down to a 
y+ = 0.5 appears necessary to obtain accurate calculations on the strength of secondary flows that 
develop in typical inlet ducts, and in particular, the size of the separated region and the engine 
face or A IP distortion signature. 

Selected use of vortex flow control to eliminate separation and improve the engine face 
distortion index can feasibly be analyzed using 3D RNS initial value procedures. The funda- 
mental nature of secondary flow, which includes internal vortex interactions that arise either from 
a vortex ingestion or vortex generators, are primarily inviscid-rotational in character, hence 
pressure driven. Thus, accurate predictions of these fundamental interactions are possible using 
simple mixing length turbulence models. 

The present investigation also demonstrates the ability of the three dimensional RNS 
techniques to analyze the flow physics associated with vortex ingestion in general geometries such 
as the F/A-18 inlet duct. Although no experminental validation data exsist on this phenomena, 
the overall correctness of the computed flow physics can be argued from simplified rotational 
inviscid flow theory. An overall computing time on the CRAY XMP of CPU = 14.5 min. for 
4.75 x 10 5 grid points makes this RNS computational approach ideally suited for routine use 
within an analysis design environment. 


18 


REFERENCES 


1. Advisory Group for Aerospace Research and Development (AGARD), "Engine Response 
to Distorted Inflow Conditions," AGARD CP-400, Sept., 1986. 

2. Bowditch, D. N. and Coltrin, R. E., "A Survey of Inlet/Engine Distortion Comparability," 
AIAA-83-1166, June 1983. 

3. Vatsa, V. N., and Wedan, B. W., "Development of an Efficient Multigrid Code for 3-D 
Navier Stokes Equations," AIAA Paper No. AIAA Paper No. 89-1791, June, 1989. 

4. Briley, W. R. and McDonald, H., "Analysis and Computation of Viscous Subsonic Pri- 
mary and Secondary Flow," AIAA Paper No. 79-1453. 

5. Briley, W. R., and McDonald, H., "Three-Dimensional Viscous Flows with Large Sec- 
ondary Velocities," Journal of Fluid Mechanics, March 1984, vol. 144, pp. 47-77. 

6. Kreskovsky, J. P., Briley, W. R. and McDonald, H., "Prediction of Laminar and Turbu- 
lent Primary and Secondary Flows in Strongly Curved Ducts," NASA CR-3388, Feb. 1981 

7. Levy, R., Briley, W. R., and McDonald, H., "Viscous Primary/Secondary Flow Analysis 
for Use with Nonorthogonal Coordinate Systems," AIAA Paper No. 83-0556, Jan. 1983. 

8. Towne, C. E., and Anderson, B. H., "Numerical Simulation of Flows in Curved Diffusers 
with Cross-Sectional Transitioning Using a Three Dimensional Viscous Analysis," AIAA 
Paper No. 81-0003, Jan. 1981. 

9. Towne, C. E., "Computation of Viscous Flow in Curved Ducts with Experimental Data," 
AIAA Paper No. 84-0531, AIAA 22th Aerospace Sciences Meeting, Jan. 1984. 

10. Anderson, B. H., "Three Dimensional Viscous Design Methodology of Supersonic Inlet 
Systems for Advanced Technology Aircraft," Numerical Methods for Engine-Airframe 
Integration, American Institute of Aeronautics and Astronautics Series, 1986. 

11. McDonald, H. and Camarata, F. J.: "An Extended Mixing Length Approach for Com- 
puting the Turbulent Boundary- Layer Development, Proceedings, Standford Conference 
of Turbulent Boundary Layers," Vol. I, Pub. by Stanford University, pp. 83-98, 1969. 

12. Reyner, T. A. and Flugge-Lotz, L, "Interaction of a Shock Wave with a Laminar 
Boundary Layer," International Journal of Non-Linear Mechanics, Vol. 3, 1968. 

13. Taylor, D., "Application of Vortex Generator Mixing Principile Principle to Diffusers," 
U.A.C. Rep. R- 15064-5, Dec. 1948. 

14. Valentine, E. F. and Carrol, R. B., "Effects of Several Arrangements of Rectangular 
Vortex Generators on the Static Pressure Rise Through a Short 2: 1 Diffuser," 

NASA RM L50L04, Feb. 1951. 


19 



REFERENCES 


15. Valentine, E. F. and Carrol, R. B., "Effects of Some Primary Variables of Rectangular 
Vortex Generators on the Static Pressure Rise Through a Short Diffuser," 

NACA RM 52B13, May 1952. 

16. Pearcy, H. H. and Stuart, C. M., "Methods of Boundary-Layer Control for Postponing 
and Alleviating Buffeting and other Effects of Shock-Induced Separation," Presented at 
the IAS National Summer Meeting, Los Angles, Calf., June, 1959. 

17. Kunik, W. G., "Application of a Computational Model for Vortex Generators in Subsonic 
Internal Flows," AIAA Paper No. 86-1458, June 1986. 

18. Squire, R., "Growth of a Vortex in Turbulent Flow," The Aeronautical Quarterly, 

Vol. 16, Pt. 3, August 1965, pp 302-306. 

19. Vakili, A. D., Wu, J. M., Liver, P. A., and Bhat, M. K., "Experimental Investigation of 
Secondary Flows in a Diffusing S-Duct with Vortex Generators," 

NASA NAG3-233, July, 1986. 

20. Horlock, J. H. and Lakshminarayana, B., "1973 Secondary Flows: Theory, Experiment, 
and Applications in Turbomachinery Aerodynamics," Ann. Rev. Fluid Mechanics, 

Vol 5, pp 247-280. 

21. Lakshminarayana, B, and Horlock, J. H., "Generalized Exprcssons for Secondary 
Vorticity Using Intrinsic Coordinates," J. Fluid Mechanics. Vol. 59, pp. 97-115. 

22. Batchelor, G. K., "An Introduction to Fluid Dynamics," Cambridge University Press, 
Cambridge, England, 1967. 


20 



Full Navier Stokes Analysis 

Figure (1) - Reduced Navier Stokes and Full Navier Stokes capabilities 
within an analysis-design engineering environment. 



Figure (2) - Paritioned geometry and grid generation for inlet ducts 
within an analysis-design engineering environment. 
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Figure (3) - Geometry definition of the Univ. Tennessee diffusing S-duct. 



Figure (4) - Computation grid for the Univ. Tennessee diffusing S-duct. 
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Figure (7) - Reduced Navier Stokes (RNS) solution showing secondary 
flow structure from vortex generator model, X/Z>,= 0.18. 



Figure (8) - Reduced Navier Stokes (RNS) solution showing limiting 
streamline signature of vortex generator region. 
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Figure (9) - Reduced Navier Stokes (RNS) solution showing limiting 

streamline signature of separation region within Univ. Tennessee 
diffusing S-duct. 



Figure (10) - Surface oil flow patterns showing separation region 
within Univ. Tennessee diffusing S-duct. 
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Figure (13) - Reduced Navier Stokes (RNS) solution showing limiting 
streamline signature of separation region. 



Figure (14) - Reduced Navier Stokes (RNS) solution showing particle 
traces associated with vortex liftoff, y /<?, = 0.03. 
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Figure (15) - Geometry definition of the F/A-18 inlet duct. 



Figure (16) - Computational grid for F/A-18 inlet duct. 
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Figure (17) - Location and geometry of vortex generators 
within F/A- 18 inlet duct. 



Figure (18) - RNS solution within F/A- 18 inlet duct showing limiting 
streamlines signature without vortex generators. 
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Figure (19) - RNS solution within F/A-18 inlet duct showing limiting 
streamlines signature with vortex generators. 


Swirl Angle 

Core Location 

Core Location 

Core Radius 

(«,) 

O’c/fi.. Zj /?,) 

o c ) 

(* uw 

0.0 

-0.551,-0.551 

0.779,225.0 

0.291 

-5.0 

-0.551,-0.551 

0.779,225.0 

0.291 

-10.0 

-0.551,-0.551 

0.779,225.0 

0.291 

-15.0 

-0.551,-0.551 

0.779,225.0 

0.291 

-20.0 

-0.551,-0.551 

0.779,225.0 

0.291 

-25.0 

-0.551,-0.551 

0.779,225.0 

0.291 

-30.0 

-0.551,-0.551 

0.779,225.0 

0.291 

-35.0 

-0.551,-0.551 

0.779,225.0 

0.291 

-40.0 

-0.551,-0.551 

0.779,225.0 

0.291 

-45.0 

-0.551,-0.551 

0.779,225.0 

0.291 


Table I - F/A-18 inlet duct vortex ingestion test cases* 


Swirl Angle 

Core Location 

Core Location 

Core Radius 

(a,) 

b’JRi, zJRI) 

(rJR„ 0 C ) 

(R„JR) 

0.0 

-0.551,-0.551 

0.779,225.0 

0.291 

10.0 

-0.551,-0.551 

0.779,225.0 

0.291 

20.0 

-0.551,-0.551 

0.779,225.0 

0.291 

25.0 

-0.551,-0.551 

0.779,225.0 

0.291 

30.0 

-0.551,-0.551 

0.779,225.0 

0.291 

40.0 

-0.551,-0.551 

0.779,225.0 

0.291 

-30.0 

-0.465, -0.465 

0.657, 225.0 

0.291 

-30.0 

-0.523, -0.523 

0.735, 225.0 

0.291 

-30.0 

-0.551, -0.551 

0.779, 225.0 

0.291 

-30.0 

-0.581, -0.581 

0.821, 225.0 

0.291 

-30.0 

-0.640, -0.640 

0.905, 225.0 

0.291 

-30.0 

-0.698, -0.698 

0.987, 225.0 

0.291 


Table II - F/A-18 inlet duct vortex ingestion test cases. 
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C- VORTEX, a, = +10.0° CC VORTEX, a, = -10.0° 

Figure (21) - RNS solution shoving vortex structure within F/A-18 
inlet duct, X/Z),= 0.0. 



C- VORTEX, a, = +10.0° CC-VORTEX, oc, = -10.0° 

Figure (22) - RNS solution showing vortex structure within F/A-18 
inlet duct, X/Z>, = 6.5. 
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Figure (23) - RNS solution within F/A-18 inlet duct showing limiting 
streamlines with ingestion of clockwise (C) vortex, 

^ = + 10 . 0 °. 



Figure (24) - RNS solution within F/A-18 inlet duct showing limiting 

streamlines with ingestion of counter clockwise (CC) vortex, 

a, = -10.0°. 
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Kinetic Energy Ratio, x 10* Kinetic Ener ?>' Rati °> x 10 







Figure (27) - RNS solution within F/A-18 inlet duct showing particle 
traces of baseline flow, y/ = 1. 14, <x I = 0.0° . 



Figure (28) - RNS solution within F/A-18 inlet duct showing particle 
traces with ingestion of counter clockwise vortex, 
y/5,= 1.14, a, = —10.0° . 
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Figure (29) - RNS solution within F/A-18 inlet duct showing particle 
traces with ingestion of counter clockwise vortex, 
y IS,~ 1.14, a, = -20.0° . 



Figure (30) - RNS solution within F/A-18 inlet duct showing particle 
traces with ingestion of counter clockwise vortex, 
y/S,= 1.14, a, = —30.0° . 
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Figure (31) - Downstream view of F/A-18 inlet duct showing 

particle traces resulting from ingestion of counter 
clockw ise vortex, y /5, = 1.14, <x s = —40.0° . 



Figure (32) - Computed trajectory of a s = —40.0° counter clockwise 
vortex relative to inlet duct centerline. 
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Figure (33) - Effect of vortex strength on impingement location 

at engine face for ingestion of counter clockwise vortex, 
rfRj = 0.551. 



Figure (34) - Effect of vortex entry point on impingement 

location at engine face for ingestion of counter clockwise 
vortex, a, = —30.0°. 


38 





NASA 

National Aeronautics and 
Space Administration 


Report Documentation Page 


1. Report No. 


NASA TM- 103703 
AIAA-91-0130 


2. Government Accession No. 


4. Title and Subtitle 

The Aerodynamic Characteristics of Vortex Ingestion 
for the F/A-18 Inlet Duct 


3. Recipient’s Catalog No. 


5. Report Date 


6. Performing Organization Code 


7. Author(s) 

Bernhard H. Anderson 


9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


8. Performing Organization Report No. 

E-5925 

10. Work Unit No. 

505-62-52 


11. Contract or Grant No. 


13. Type of Report and Period Covered 
Technical Memorandum 


14. Sponsoring Agency Code 


15. Supplementary Notes 

Prepared for the 29th Aerspace Sciences Meeting sponsored by the American Institute of Aeronautics and 
Astronautics, Reno, Nevada, January 7-10, 1991. 


16. Abstract 

u *- . 

A Reduced Navier Stokes (RNS) solution technique has been successfully combined with the concept of partitioned 
geometry and mesh generation to form a very efficient 3D RNS code aimed at the analysis-design engineering 
environment. Partitioned geometry and mesh generation is a pre-processor to augment existing geometry and grid 
generation programs which allows the solver to (1) recluster and existing gridfile mesh lattice, and (2) perturb an 
existing gridfile definition to alter the cross-sectional shape and inlet duct centerline distribution without returning 
to the external geometry and grid generator. The present results provide a quantitative validation of the initial 
value space marching 3D RNS procedure and demonstrates accurate predictions of the engine face flow field, 
with a separation present in the inlet duct as well as when vortex generators are installed to supress flow 
separation. The present results also demonstrate the ability of the 3D RNS procedure to analyze the flow physics 
associated with vortex ingestion in general geometry ducts such as the F/A-18 inlet. At the conditions 
investigated, these interactions are basically VMnviscid like**, i.e.jhe dominant aerodynamic characteristics have 
there origin in inviscid flow theory. 


17. Key Words (Suggested by Author(s)) 

Applied aerodynamics 
Computational fluid dynamics 
Internal flow 


18. Distribution Statement 

Unclassified — Unlimited 
Subject Category 20 


19. Security Classif. (of this report) 
Unclassified 


20. Security Classif. (of this page) 
Unclassified 


21 . No. of pages 

39 


22. Price* 


NASA FORM 1626 OCT 66 


*For sale by the National Technical Information Service, Springfield, Virginia 22161 





